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f— ^ I We present a direct measurement of the parity-violation parameter Ac in 

^D I the coupling of the Z^ to c-quarks with the SLD detector. The measurement 

D I is based on a sample of 530k hadronic Z^ decays, produced with a mean 

Q I electron-beam polarization of |Pe| = 73%. The tagging of c-quark events 



is performed using two methods: the exclusive reconstruction of D*~^, D~^, 



kS ' and D^ mesons, and the soft-pions (vr^) produced in the decay of D*^ 



D^nf . The large background from D mesons produced in B hadron decays is 
separated efficiently from the signal using precision vertex information. The 
combination of these two methods yields Ac = 0.688 ± 0.041. 
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I. INTRODUCTION 

In the Standard Model, the Z^ couphng to fermions has both vector [v^] and axial- 
vector iaf) components. Measurements of fermion asymmetries at the Z° resonance probe 
a combination of these components given by 

A^ = 2vja^l{v)^a]). (1) 

The parameter Aj expresses the extent of parity violation at the Zjj vertex and its mea- 
surement provides a sensitive test of the Standard Model. 

At the Born-level, the differential cross section for the reaction e^e~ ^ Z^ —>■ f f is 

af{z) = daf/dz oc (1 - A^P^){1 + z^) + 2A^{A^ - Pe)z, (2) 

where Pe is the longitudinal polarization of the electron beam (Pg > for net right-handed 
polarization) and z = cos 6, 9 being the polar angle of the outgoing fermion relative to the 
incident electron. In the absence of electron beam polarization, the parameter Af can be 
extracted by isolating the term linear in z via the forward-backward asymmetry: 

a} ^^{^) ~ <^^{—z) 2z 

which also depends on the initial state electron parity- violation parameter A^.. At the SLAG 
Linear Collider (SLC), the ability to manipulate the longitudinal polarization of the electron 
beam allows the isolation of the parameter Af in Eq. |^, independent of Ag, using the left-right 
forward-backward asymmetry: 

Aijz) = ^4^^ T^ -^^ T^ = PeAf -, (4) 

^^^^ [ai{z)+ai{-z)] + [aj,{z)+ai{-z)] '^'^1 + .^' ^^ 

where indices L, R refers to Z^ —>■ ff decays produced with left-handed or right-handed 
polarization of the electron beam, respectively. For a highly polarized electron beam with 
|Pe|=73%, ApQ provides a statistical advantage of (Pe/Ag)^ ~ 24 in the sensitivity to Af 
relative to the unpolarized asymmetry. 



In this paper, we present a direct measurement of the parity-violation parameter Ac for 
the Zee couphng. The c-quarkQ is the only up-type quark which can be identified, and its 
measurements provides sensitive test of the standard model. The tagging of c-quarks is 
performed using exclusively reconstructed D*'^, D~^, and D^ mesons, as well as an inclusive 
sample of D*^ -^ D'^T^t decays identified by the soft-pion (tTs). 

The charge of the primary c-quark is determined by the charge of the D^*\ K (in the 
D^ reconstruction case), or vr^ (in the soft-pion analysis case). The direction of the primary 
quark is estimated from the direction of the exclusively reconstructed /)(*)+ or D^ meson, or 
the jet axis in the soft-pion analysis. The value of A^. is extracted via an unbinned maximum 
likelihood fit. The fit is performed on two separate data samples: one collected between 1993 
and 1995, and the other, with an improved vertex detector, between 1996 and 1998. The 
data samples associated with these two periods comprise 150k and 380k hadronic Z° decays, 
respectively. 

The measurement of Ac presented in this paper updates and supersedes our previous 
publication [l[], which was based on a sample of 50k hadronic Z^ decays from 1993 alone. 
There are several direct and indirect Ac measurements ||T]-§1. The measurement reported 
here is currently the most precise. 

II. APPARATUS AND EVENT SELECTION 

The measurement described here is based on 530k hadronic Z^ decays recorded in 1993- 
98 with the SLC Large Detector (SLD) at the SLC e^e~ collider at a mean center-of-mass- 
energy of 91.27 GeV(1993-95) or 91.24 GeV(1996-98). A general description of the SLD can 
be found elsewhere [Q. Charged-particle tracking for the 1993-95 data sample uses the cen- 
tral drift chamber (CDC) and VXD2 CCD pixel vertex detector. For this system, the 
measured impact-parameter resolution in the transverse (longitudinal) direction with respect 



^Throughout the paper charge-conjugate states are implicitly included. 



to the beam axis can be approximated by 11 70/Psin^'^ 6 /im (38 70/Psin^'^ 9 /xm), as 
a function of the track momentum P (in GeV/c) and the polar angle 6. In 1996, we installed 
the upgraded 307M pixel vertex detector (VXDS) ||^, which provides improved impact- 
parameter resolution of 7.8 33/P sin^' ^ 9 /im (9.7 33/P sin^' ^ 9 yum) [^ in the transverse 
(longitudinal) direction with respect to the beam axis. In addition, VXD3 extended the 
polar-angle coverage from |cos0| < 0.75 to |cos0| < 0.85. Combining the CDC and the 



VXD, a momentum resolution of a{PT)/PT = v (-01)^ + {.QQ2QPt / GeVy is achieved. The 
Liquid Argon Calorimeter (LAC) [Q measures the energy of charged and neutral particles 
and is also used for electron identification. The barrel LAC covers the polar-angle region of 



I cos6'| < 0.84, and has energy resolutions of 15%/wP(GeV) and 65%/wP(GeV) for electro- 
magnetic and hadronic showers, respectively. Muon identification is provided by the Warm 
Iron Calorimeter (WIC) [|TU|. The Cerenkov Ring Imaging Detector (CRID) |[TT]] provides 
particle identification. In order to achieve particle identification over a wide momentum 
range, the CRID uses two different radiator systems; liquid (C6F14) and gas (C5F12), which 
provide excellent tt / K separation in the momentum range from 0.3 to 35 GeV/c. 



The SLC operates a polarized electron beam and an unpolarized positron beam ||12 



The average electron polarization measured for the 1993-98 data sample is |Pe| = 73 ± 0.5% 



12,1^. The SLC interaction-point (IP) size in the xy plane is 2.6 /im x 0.8 /xm and its mean 



position is reconstructed with a precision of aip = 4/im (7/im) using the tracks in sets of 



~30 sequential hadronic events for the 1996-98 (1993-95) data sets [14|. The event-by-event 
median z position of tracks at their point of closest approach to the IP in the xy plane 
determines the z position of the Z^ primary vertex (PV) with a precision of ~15 /im (35 
/im) for the 1996-98 (1993-95) data. 

Hadronic events are selected by requiring at least 5 charged tracks, a total charged energy 
of at least 20 GeV/c, and a thrust axis calculated from charged tracks satisfying | cos 9thrust\ < 
0.87 (0.8 for the 1993-95 data). In the event selection and charm reconstruction, we use the 
quality tracks which satisfy the following criteria for the 1996-98 (1993-95) samples: 



1. At least 23 (30 for the 1993-95 data) associated CDC hits; 

2. A radius of the innermost CDC hit of the reconstructed track within 50 cm (39 cm) of 
the IP; 

3. An xy and rz impact parameter with respect to the IP of less than 5cm (10cm); 

4. A reconstructed polar angle within | cos 6* |< 0.87 (0.80); and 

5. A momentum component transverse to the beam axis greater than 0.15 GeV/c. 

As Z^ -^ bb events are also a copious source of D mesons, they represent a potential 
background. We reject these events using the invariant mass of the charged tracks associ- 



ated with the reconstructed secondary decay vertices [T^. In particular, we require that 
there must be no vertex with invariant mass greater than 2.0 GeV/c^. Monte Carlo (MC) 
simulations indicate that this cut rejects 57% of bb events while preserving 99% of cc events. 

III. Ac MEASUREMENT WITH EXCLUSIVE CHARMED-MESON 

RECONSTRUCTION 

In this analysis, we reconstruct three different D^*^ meson states for c-quark tagging: 
the pseudo-scalar mesons D^ and D^, and the vector meson D*^. This section describes 
the procedure for their reconstruction, as well as the corresponding Ac measurement and a 
discussion of associated systematic errors. 

A. Z)*+ selection 
D*^ mesons are identified via the decay D*^ -^ D^T^t followed by: 

D° -^ K-n+TT^ Satellite, 

D^ -^ K^7i^7i~7i~^ K-n-RTi, or 

D^ -^ K^l^ui (/=e or /x) Semileptonic. 



In these decays, the charge of the underlying c quark is specified by the charge of the "soft 
pion" TTg. No attempt is made to reconstruct the tc^ in the sateUite mode, nor to estimate 
the neutrino direction or energy in the semileptonic mode. 

We search for D*^ mesons in each of the two event hemispheres, defined by the plane 
perpendicular to the thrust axis, using all quality tracks with at least one hit in the VXD. 
In the A'tttttt mode, we only use tracks which have momentum greater than 0.75 GeV/c. We 
first construct D^ candidates using all combinations of tracks corresponding to the charged 
multiplicity in each Z}° decay mode, with zero net charge. Here one of them is assigned the 
charged kaon mass and the other(s) are assigned the charged pion mass. In the semileptonic 
mode, we combine an identified electron or muon track with another track which has opposite 
charge and assume the track to be a kaon. Electrons are identified based on the momentum 



measured with the CDC and the energy deposited in the calorimeter [16|. Electrons from 



7-conversions are rejected. Muon candidates are identified by the association of extrapolated 
CDC tracks with hits in the WIC [0]. 

A vertex fit is performed on the tracks forming a D^ candidate, and we require that its 
X^ probability be greater than 1%. The invariant mass M of the D^ candidates is required 
to lie within the following ranges: 

1.765GeV/c^ < M^o < 1.965GeV/c^ {Ktt), 
1.500GeV/c^ < Mdo < 1.600GeV/c^ (Satellite), 
1.795GeV/c^ < M^o < 1.935GeV/c^ (KTnnr), 
l.lOOGeV/c^ < Mdo < 1.800GeV/c^ (Semileptonic). 

These reconstructed pseudo-scalar meson candidates are then combined with a soft-pion 
candidate track with charge opposite to that of the kaon candidate, thus forming the D*~^ 
candidate. 

To reconstruct the D*^, we use two sets of selection criteria. One is based on event 
kinematics and the other on event topology. The former relies on the fact that D*^ mesons 
in cc events have much higher xd* = '2E£)*/Ecm, where E^* is the D*~^ energy, than those 



in bb events or random combinatoric background (RCBG). The latter relies on the fact that 
D°'s in cc events have a longer 3D decay length (~ 1 mm) than that for RCBG, and originate 
at the primary vertex, in contrast to those Z^^'s in bb events emerging from B decay vertices. 
We select the combinations which satisfy either condition. 

In the selection based on the event kinematics, we require the candidate to have xd* 
greater than 0.4 {Ktt, Satellite, and Semileptonic) or 0.6 {Ktttiti). For a true D^ candidate, 
the distribution of cos^*, where 6* is the opening angle between the direction of the D^ in 
the laboratory frame and the kaon in the D^ rest frame, is expected to be flat. Since 
background events peak at cos^* = ±1, they are further reduced by requiring | cos^*| < 0.9 
[Kir, Satellite, and Semileptonic) or 0.8 {Ktttiti). We also require the soft-pion candidate to 
have momentum greater than 1 GeV/c. In the satellite mode, we apply a 3D decay-length 
cut of L/ai > 1.5 on the reconstructed D^ vertices to reduce the RCBG. (The average 
decay-length resolution is < (Tl > ~ 150 /im.) 

In the selection based on the event topologies, we require the reconstructed D^ vertices 
to have 3D decay-length significance L/ctl > 2.5, and the xy impact parameter of the D^ 
momentum vector to the IP to be less than 20 /xm {Ktt and Ktttiti) or 30 fira {Satellite and 
Semileptonic). The latter cut is effective in rejecting D decays in bb events. Since these D's 
have significant Pt relative to the parent B flight direction, and the 5's themselves have 
a significant flight length (~ 3.5 mm), many of these D's do not appear to originate from 
the primary vertex. A cut of xd* greater than 0.3 {Ktt, Satellite, and Semileptonic) or 
0.4 {Ktttttt) is also applied. Fig. 1 shows the distribution of xy impact parameter of the D^ 
relative to the IP for the decay of D*^ -^ D^tt^ , D^ -^ K^tt^ . In this flgure, we do not 
reject -B-decay candidate events with the invariant mass cut of the reconstructed secondary 
vertices described above, only for the purpose of showing how the xy-impact-parameter cut 
is effective in rejecting the -B-decay background. After applying the invariant mass cut of 
the reconstructed secondary vertices, 34% of the remaining 5-decay background events are 
rejected by the xy-impact-parameter cut. 



The overlaps of the sets of candidates from the event kinematics and topology analysis 
are 53% (Kn), 50% (satellite), 28% {KnTCTr), and 36% (semileptonic). In the Kninr sample, 
there may be multiple D^ candidates in a single event which pass the above cuts. To avoid 
double counting and to reduce the background, we select the D^ candidate with the lowest 
vertex x^- 

Having selected a candidate, we form the mass difference AM = Md* — M^o. The 
mass difference spectra for the four reconstructed D*^ decay modes are shown in Fig. 2. 
For all decay modes, clear peaks around AM = 0.14 GeV/c^ appear due to the D*^ to 
D^ transition. We include the candidates in the signal sample provided AM is less than 
0.148 GeV/c^ {Kn and KTrmr), less than 0.155 GeV/c^ (Satellite), and less than 0.16 
GeV/c^ {Semileptonic) . The side-band region is defined as 0.16 < AM < 0.20 GeV/c^ 
(0.17 < AM < 0.20 GeV/c^ for the Semileptonic mode), and is used to estimate the RCBG 
contamination in the signal region. In the figure, the MC predictions for the reconstructed 
D*^ (open) and RCBG (hatched) are also presented. For the MC prediction, the relative 
normalizations of signal and RCBG shapes are adjusted so that the predicted numbers of 
events match those observed in the data signal and side-band regions. Averaged over the 
various modes, this procedure requires adding 10% to the MC signal and 5% to the MC 
RCBG. The number of the selected candidates as well as the contributions oi c,b ^ D and 
RCBG estimated by MC are summarized in Table |l[ 

B. D^ and D^ selection 

The D~^ and D^ mesons are identified via the decay channels 

D+ -^ K-TT + TT + 



D'^^K- 



n ' 



These modes are reconstructed by considering all quality tracks in each hemisphere which 
have VXD hits. In the D~^ reconstruction, we additionally require each track to have a 
momentum of greater than 1 GeV/c. 



For the D^ reconstruction, we combine two same-sign tracks, assumed to be pions, with 
an opposite-sign track, assumed to be a kaon. We require that xd+ be greater than 0.4, and 
cos 6'* be greater than -0.8, where 9* is the opening angle between the direction of the D^ in 
the laboratory frame and the kaon in the D^ rest frame. To reject D*^ decays, the differences 
between Mk--k+-k+ and Mk--k+ are formed for each of the pions, and both are required to be 
greater than 0.16 GeV/c^. To remove RCBG, we require that the x^ probability of the good 
vertex fit be greater than 1%, and that the 3D decay-length significance L/a^ be greater 
than 3.0. To reject D~*"'s from hb events, the angle between the D^ momentum vector and 
the vertex flight direction is required to be less than 5 mrad in xy and less than 20 mrad 
in rz. Here we use the angular information instead of the impact-parameter information. 
We can strongly constrain the D^ to originate from the IP with the angular information, 
because of its large decay length. 

To form the D^ vertices, tracks identified as charged kaons, by the requirement that the 



GRID log- likelihood [17] for the K hypothesis exceeds that for the vr hypothesis by at least 3 
units, are combined with an opposite-charge track, assumed to be a pion. We use the GRID 
information for this mode only. To reject background we require x^o be greater than 0.4. 
We require that the vertex fit have x^ probability greater than 1% and the 3D decay-length 
cut L/ai be greater than 3.0. To reject the D°'s from D*^ decays, the differences between 
Mk--k+-k+ or MK-Tt+Tz--, and Mk--k+ are formed for all other tracks in the same hemisphere, 
and these are required to be greater than 0.16 GeV/c^. Finally, to reject Z^^'s from hh events, 
we require that the xy impact parameter of the Z)" momentum vector relative to the IP be 
less than 20 /im. 

D+ and D° candidates in the ranges of 1.800 < MK--n+iT+ < 1-940 GeV/c^ and 1.765 
< Mk-tt+ < 1.965 GeV/c^, respectively, are regarded as signal. The side-band regions are 
defined as 1.640 < Mk~^+^+ < 1.740 GeV/c^ and 2.000 < Mk-^+^+ < 2.100 GeV/c^ for 
D+, and 2.100 < Mx-t,+ < 2.500 GeV/c^ for D°. In Fig. 3, the invariant mass spectra 
for the resulting D^ and D^ signals are plotted. The backgrounds in the signal regions are 



estimated from the MC in the same manner as in the D*^ analysis. 

C. Measurement of Ac 

Using the six decay modes, we select 3967 D*+, /}+, and D^ candidates from 1993-98 SLD 
data. The estimated composition is 2829±35 c —>■ D signal, 281±11 b ^ D, and 857±19 
RCBG. These c ^ D signals correspond to a selection efficiency for cc events of 3.9%. The 
results for the number of selected candidates are summarized in Table J[ 

The charge of the primary c-quark is determined by the charge of the D^*\ or K (in 
the D^ case). The direction of the primary quark is estimated from the direction of the 
reconstructed D meson. Fig. 4 shows q cos Od distributions, for the selected D meson sample 
separately for left- and right-handed electron beams. Here, q is the sign of the charge of the 
primary c-quark and 6d is the polar angle of the reconstructed D meson. 

To extract Ac, we use an unbinned maximum likelihood fit based on the Born-level cross 
section for fermion production in Z'^-boson decay. The likelihood function used in this 
analysis is 



\nC = Y: ln{Pi(a;i,) ■ [(1 - PcAc){l+y^) + 2{Ac - Pe)yi ■ Af] 
1=1 

+ Pl{x'^) ■ [(1 - PcAc){l + yl) + 2{Ac - Pe)y.. ■ Af ] 

+ PkcBci^D) ■ [(1 + yl) + 2A^CBGl/d} (5) 

where y = qcosOr,, n is the total number of candidates, and the index j indicates each of the 
six charm decay modes. 

A^ and A^ are the asymmetries from D*+, D^ , and D^ mesons in cc and hb events, 
respectively. We treat A^ as a free parameter, while A^ is fixed. A^ is estimated in a 
similar manner to Ref. jT^. We start with the Standard Model prediction |]TP|, Af, = 0.935, 
and assign it an error of ±0.025 from the average value of SLD measurements of 0.911 ± 
0.025 PD|. This 6-quark asymmetry is diluted by B^-B^ mixing and the wrong-sign D meson 



from the W in b —>■ cW , W —>■ cs decay. The effective b asymmetry can be expressed by 
correcting with two dilution factors: 

A^ = AbX (1- 2xmixing)C^ - 2xw-^cs)- (6) 

The value of Xmixing is deduced from the D-meson production rates through B decays. We 
estimate the B ^ D source fractions from MC Using the fractions and the x values of 
X = 0.1186 ± 0.0043 and Xd = 0.156 ± 0.024 [^, we derive the Xmixing value for D*+, 



D^, or D^. The value of Xw-->cs, the correction for wrong-sign D mesons from the W~ 
in 6 — » cW~ decay, is also estimated from MC. We obtain xw-^cs = 0.023 ± 0.006 for the 
average of D*~^, D~^, and D° mesons, and 0.021 ±0.006 for D*~^ mesons only. Here the errors 
include the theoretical error of 30% coming from Br{b — ^ ccs) = 22 ± 6% p2|. The former 
and latter xw-^cs values are used for exclusive D reconstruction and inclusive soft-pion 
analysis, respectively. By combining these two dilutions, we obtain 

A^ = 0.657 ±0.025 for £)*+, 
= 0.655 ± 0.026 for D+ and, 
= 0.762 ± 0.023 for D^. 

To check the Af value, we measure Aj^ for D*~^ using the 1996-98 experimental data. 
In this measurement, we select D*~^ mesons in the decay, D*^ -^ D^^^t followed by D^ — >■ 
K^-K^ , D^ -^ K^TT^TT^, or D^ -^ K^ti^ti^ii^ . The bb events are selected by requiring that 
the invariant mass for the reconstructed secondary vertices be greater than 2 GeV/c^ for 
at least one of the two event hemispheres. In order to select the D*~^ mesons, we apply 
similar cuts to those used to select the D*~^ mesons from c-quarks, but without any xy 
impact parameter cut to reject D*~^^s from 6-quarks. We select 2196 D* candidates with 
the fractions of 63% b ^ D, 2% c ^ D, and 35% RCBG. Using this sample, we measure 
Af = 0.58 ± 0.10, which is consistent with our assumed Aj^ value for D*^. The error of 0.10 
is treated as a systematic error of Aj^. 

We also check the effect of the decay-length cut of the reconstructed D mesons. In this 
analysis, we apply the decay-length cut of L/a^ > 1.5 ~ 3.0 (depending on the charm decay 



mode) to reject RCBG. This cut may increase the effective value of Xmixing- Using our MC, 
we estimate the effect of this cut to be small {Axmixing / Xmixing = 3%). 

Arcbg is the analog of Ac for the RCBG, and we expect it to be very small. The 
asymmetry in the side-band region is measured as —0.0006 ± 0.0031, and is assumed to be 
zero. For Ae, we have taken A^ = 0.1513 ± 0.0022 from the SLD measurement | ]13|| . 

P/, Pj^, and Prcbg ^^^ ^^^ probabilities that a candidate from the jth decay mode is a 
signal from cc, bb, or RCBG. The determination of these functions is based on the relative 
fractions and the xd distributions for the six decay modes. They are defined as: 

Nsignalixo) fcixo) 



Pc{Xd) 



NtotaliXo) fciXo) + fbiXo) 



p , X _ NsignaliXp) _ fbjXp) . . 

Ntotaiixn) fcixo) + fbixn) 

where Ntotaiixo) is the observed number of D mesons, and Nbg{xd) is that of background 
events, in the x^ bin. Using the xjj distributions for the reconstructed D mesons and side- 
band events, we determine the ratio NBc/Ntotai in each xd bin. The ratio N signal /N total is 
given by the relation N signal /N total = 1 - NsclNtotai in each bin. Figs. 5 (a)-(f) show the xd 
distributions for six decay modes, which are used in this determination. 

The functions fcixo) and fbixo) describe the fraction of D mesons in the c and b decays, 
respectively, and are expressed as 

fc{b) = Wc(fe) ■ dc(b)ixD), (8) 

where dc{b){xD) describes the shape of xc distributions in c(6) -^ D, and ujc(b) represents 
the total fraction of the c{b) -^ D for the reconstructed D candidates. We obtain the 
function dc{b){xD) from MC, and the values of Uc and ujfy are derived from Table |. The ratio 
fc{b)/{fc + fb) gives the probability that a D candidate is from a primary c{b) quark. 

Performing the maximum likelihood fit to the data sample, we measure Ac = 0.671 ±0.096 
(1993-95) and Ac = 0.681 dz 0.047 (1996-98). As a check, we also determine Ac with a simple 



binned fit of tlie type described in Ref. ^. We find A^ = 0.731 ± 0.102 (1993-95) and 



Ac = 0.666 ± 0.049 (1996-98); which are consistent with the values above. 

D. QCD and QED correction 

As a result of hard gluon radiation, the extracted value of Ac(b) is somewhat different than 
its Born-level value in Eq. |l|. To account for this, the fit parameter ^^(fc) in the likelihood 
function is replaced with the first-order corrected parameter Ac(b)i^ — ^qcd(cos6')) with 
AqJ^(cos^) = Cc(t,)AgJ^ 5o(cos6'), where AgJ^ indicates the magnitude of the leading- 
order (LO) QCD correction for c(6)-quark production, and ^qcd so i^ ^^^ LO QCD cor- 



rection calculated by Stav and Olsen including the quark-mass effect [^. The factor Cc(b) 
takes into account the mitigation of the effects of gluon radiation due to the analysis pro- 
cedure. For example, the requirement that D mesons have high xd values selects against 
events containing hard gluon radiation, reducing the overall effect of gluon radiation on the 
observed asymmetry. 

The correction factor Cc(b) is estimated with the MC by comparing the effects of QCD 
radiation, for the JETSET Parton Shower model, with and without the full analysis including 
detector simulation: 

/Iffen Ameas 

^g ~ A gen A gen [Q ~ '^y'^Jy V^) 

where the superscripts '(/en' and ^meas^ refer to the MC asymmetries for generator level 
(Parton Shower Model simulation only) and fully analyzed events, respectively. These MC 
asymmetries are determined by doing a fit to the form 

, 2 cos 9 , , 

in bins of cos^. We obtain Cc = 0.27 ± 0.10 and Cb = 0.17 ± 0.08 for c-quark and 6-quark, 
respectively. Applying the first-order QCD correction with the correction factors Cc{b), leads 
to a 1.0% increase of Ac. 



In this analysis, we have also considered the effects of next-to-leading order (NLO) gluon 
radiation. The NLO QCD correction is written as: 

2 



A?^"'^= - x4.4xCe + A,,. (11) 



where the first term is from hard gluon emission [|J]. We use the same correction factor C^ 
as in Eq. |. The second term A^s accounts for the effects of the process g —^ cc for gluons 
which arise during the shower and fragmentation processes. 

The effects of gluon splitting have been taken into account by analyzing the MC as if 
it were data, with and without events with gluon splitting. The resulting difference must 
be scaled to account for the difference between the JETSET gluon splitting rates and the 
currently measured values for these rates. The rate for gluon splitting to charm quark pairs 
in JETSET is 0.0136 per hadronic event, and the current LEP average |2^ is 0.0319 ± 
0.0046, yielding a scale factor of 2.35 ± 0.34. 

The second-order QCD correction increases Ac by 0.4%. Applying the first- and second- 
order QCD corrections, we obtain A^ = 0.681 ± 0.097 (1993-95) and A^ = 0.690 ± 0.047 
(1996-98). 

Using ZFITTER(6.23) |T9|, we estimate QED corrections including initial- and final- 
state radiation, vertex correction, 7 exchange, and 7-Z interference. We use the input 
values nitop = 175 GeV/c^ and mmggs = 150 GeV/c^. These corrections increase Ac by 0.2%. 
Applying the QED corrections, we obtain Ac = 0.682±0.097 (1993-95) and Ac = 0.691±0.047 
(1996-98). 

E. Systematic errors 

The following systematic errors have been estimated and are summarized in Table |I|: 

• The largest uncertainties are due to the RCBG, arising from the statistics of the MC 
and side-band events, which are used to determine the fraction of the RCBG in the 
signal, and the shape of RCBG X£, distribution which is determined by side-band 



events. The uncertainty of the RCBG xd shape is estimated by comparing the xd 
distributions for MC RCBG events and for side-band events. 

There is a difference in acceptance between signal and RCBG event samples. In this 
analysis, we determine the RCBG probability function as a function of xd- This 
is correct if the ratio between the signal and RCBG acceptance is constant over the 
different cos 6 regions. In order to study this, we compare the RCBG | cos ^| distribution 
obtained from the side-band region and that from the signal region events weighted by 
the RCBG probability function Prcbg{^d) in Eq- ^ These two distributions become 
significantly different starting at |cos0| ~ 0.65. Hence, we apply an acceptance cut 
of I cos^dI < 0.65, then regard the difference between with and without the cut as a 
systematic uncertainty. 

We expect the asymmetry of RCBG to be very small, and take a central value of 
Arcbg = 0. Since the asymmetry of the side-band events is measured to be — 0.0006 ± 
0.0031, we take -0.0037 as a lower limit on Arcbg- 

We vary fb^n/ifb^D + fc^n), the fraction of D mesons from Z° -^ bb, by ±20% 
to account for differences between our MC and the range of measurements of /}(*)+ 
production in Z° decay [ITSJ ^ 



The effect of the uncertainty of A^ is estimated by varying 6Aj^ = ±0.10, where the 
error is from the statistical error of our Aj^ measurement by using experimental data. 
In Table |I|, we show the resultant error in A^ coming from the uncertainty in Af, 
(0.935 ± 0.025) separately from the uncertainty in the mixing parameter. 

The systematic error on the fragmentation function is estimated by modifying the 
Xd distributions in heavy-quark fragmentation. In our MC sample, we use Peterson 
fragmentation and the average xd values are (xd) = 0.508 and 0.318 for c ^ -D and 
b —>■ D, respectively. We change the values by A{xd) = ±0.015(±0.010) for c{b) -^ D. 



Our sensitivity to the RCBG xd distribution is checked by performing the analysis 
with Prcbg derived from the MC background instead of the data side-bands. 

The shapes of the xd distributions in c{b) -^ D, expressed as dc(b)ixD) in Eq. H, are 
obtained by fitting to the MC xo distributions. The sensitivity to this procedure is 
checked by performing the analysis with a binned MC xd distribution. 

We assume Ag = 0.1513 ± 0.0022, and estimate this systematic error by varying Ae 
within the error. The precision of the polarization measurements are APg = 1-1% 
(1993), 0.5% (1994-95), and 0.4%(1996-98) 0,0]. We estimate the systematic error 
due to polarization uncertainties by varying Pg with these errors. 

We consider two sources of uncertainties on the leading order QCD correction: The 
uncertainty on a^ and the uncertainty in the estimation of the correction factor due 
to the analysis bias. The range of a^ chosen for the analysis is 0.118 ± 0.007, while 
that for the correction factor is 0.27 ± 0.10 for c-quark or 0.17 ± 0.08 for 6-quark, as 



described in Section HID. 



• In order to estimate the hard-gluon-radiation uncertainty in the second-order QCD 
correction, we vary the magnitude of the correction by 50% of itself. We use the 
experimental error for the uncertainty in gluon splitting into cc. 

The total systematic errors are 0.034 and 0.021 for 1993-95 and 1996-98 SLD runs, re- 
spectively. 

F. Results 

We obtain the following results for the measurements using exclusive channels: Ac = 
0.682 ± 0.097(stat.) ± 0.034(s7/s.) (1993-95) and A = 0.691 ± 0.047(stat.) ± 0.02l{sys.) 
(1996-98). The combined result is: 

Ac = 0.690 ± 0.042(stot.) ± 0.021{sys.) 



IV. INCLUSIVE SOFT-PION ANALYSIS 

In this analysis, c-quarks are identified by the presence of soft pions from the decay 
D*^ -^ D^Ti^ . Since this decay has a small Q value of mr,* — m£,o — mj^ = 6 MeV/c^, the 
maximum transverse momentum of the vr^ with respect to the D*^ flight direction is only 40 

MeV/c. 

A. Jet reconstruction and soft-pion selection 

We select hadron events and reject bb events by using the same criteria described in 
Section 0. The D*'^ flight direction is approximated by the jet direction, where charged 
tracks and neutral clusters are clustered into jets, using an invariant- mass (JADE) algorithm. 
In the jet clustering, particles are merged together in an iterative way if their invariant mass 
is less than 4.6 GeV/c^. We only use the tracks and clusters which have the momentum of 
greater than 1.2 GeV/c and 1.0 GeV/c, respectively, to form the jet. The tracks are required 
to satisfy the track quality cuts described in Section || and to have vertex hits. 

The jets must satisfy the following criteria: 

1. At least 3 charged tracks; 

2. At least one track with momentum P > 5 GeV/c; 

3. The net charge of the jet. Eg, should be |Sg| < 2; 

4. Sum of the largest and second largest 3D normalized impact parameters of the tracks 
> 2.5 cr; and 

5. There is at least one opposite-charged-track pair which has x^ probability of two tracks 
coming from the same vertex greater than 1%. 

The criteria 2) and 3) are effective to reduce the huge RCBG. The criterion 4) rejects the 
light flavor events. The criterion 5) relies on the fact that it is likely that the D^ decays into 
at least one pair of oppositely charged tracks. 



After selecting the jet candidates, we look for the soft-pions using a momentum cut of 
1 < P < 3 GeV/c and an impact-parameter cut of less than 2a from the IP. Since soft-pions 
in cc events have much higher momentum than those in hh events, the former criterion rejects 
such soft-pions from hb events. The latter criterion is also effective to reduce the soft-pions 
from 66, because D* decays from hh events have significant transverse momentum relative to 
the parent B flight direction, and they do not appear to originate from the primary vertex 
due to the B lifetime. 

Using the selected soft-pion candidates, the momenta transverse to the jet axis. Pt-, are 
calculated. Fig. 6 (a) shows the P^ distribution for the soft-pion candidate tracks. The peak 
around P^ = is from charm signal. We define P^ < 0.01 (GeV/c)^ as the signal region, 
where a signal-to-background ratio of 1:2 is observed. From 1993-98 data, 12992 soft-pion 
candidates are selected in the region. 

B. BG determination and Ac measurement 

To evaluate the number of the D*^ -^ D^^^t decays, a fit to the observed P^ distribution 
is performed using the signal plus background shape. The signal shape is assumed to be a 
simple exponential 

S{P^) = aexp{-P^/(3). 

We obtain (3 = 0.00471 ± 0.00007 by fitting the MC spectrum of D*+ -^ P)°7r+ decays and 
fix the value of f3 to fit the experimental data. For the background shape, we try two kinds 
of functions with three free parameters each: 

Pi(P|) = a/(l + 6P| + c(P|)'), 
P2(p2) = a' + 6'exp(-P|/c'). 

The fit results are illustrated in Fig. 6 (a), where we show the extrapolation of Pi(P^) 
(dashed line) and F2{P^) (dotted line). 



The observed signal in 1993-98 data is 4291 ± 147 (xV^df = 219.0/196) with 5(P|) + 
Fi(P|) and 4032 ± 124 (xVndf = 224.0/196) with 5(P|)+F2(P|), where the fit is performed 
in each case for P^ < 0.1 GeV/c. We choose Pi(Py) for the background shape to measure 
the Ac, because of its smaller x^/ndf value. The difference between these two functions is 
regarded as a systematic error. 

We determine the relative normalizations of signal and background for the MC prediction 
using the above fit to the data. Fig. 6 (b) shows the detailed P^ distribution from the MC 
prediction with this normalization. We also overlay the background shape extrapolated by 
the fitting with S{P^) +Pi(P|) (dashed line). Using the MC, we estimate the contributions 
of c -> D*+ and b -^ D*+ as 3791±39 and 500±14, respectively, in 1993-98 data. 

In order to ensure that there is little room for non-D* sources of slow pions in the data, 
we compared the signal obtained by fitting to the experimental data and the number of D*'s 
expected by MC. Here normalization of the MC is determined by the number of hadronic 
events. Using MC, we estimate the number to be 4507±57. Comparing this number and the 
obtained experimental number of 4291 ± 147, we conclude that other charm-decay sources 
in the experimental data are small. 

The direction of the primary quark is estimated from the jet axis, and the charge of the 
primary c-quark is determined by the charge of the vr^. Fig. 7 shows the q cos Od distributions, 
where q is the sign of the primary c-quark, and Oo is the polar angle of the jet axis, for the 
selected D*^ sample separately for left- and right-handed electron beams. 

To extract Ac-, we use an unbinned maximum likelihood fit, using a likelihood function 
similar to the exclusive D reconstruction analysis (Eq. ^). We regard the Ac as a free 
parameter, and fix the asymmetry of D*'^ from hh events, A^ . This value is obtained by 
following the similar procedure described in Section [1IIC|) . 

We expect the asymmetry for the BG, Abg-, to be very small and assume it to be zero. 
Using the MC, we measure the asymmetry of the background to be 0.009 ± 0.017. 

For the probabilities Pc, Pb, and Prcbg in Eq- ^ we used the following functions: 



p/p p2\ _ ^^signaiy-T, rrp) Jc\^) 



Ntotal{P,P^) UP) + h{P) 
Pb{P,PT) 



,2\ Nsignal{P, Pt) fb{P) 



Ntotal{P,P^) fc{P) + MP) 
p (p p2N _ Nbg{P,Pt) . (-,^^ 



where P and P^ indicate the momentum and the squared transverse momentum to the D* 
jet axis for soft-pion tracks, respectively. Ntotai and Nbg are the observed number of soft- 
pion candidates and that of background in each P and P^ bin, respectively. We estimate 
Nbg from MC, and the relation N signal /Ntotai = ^- Nbg /Ntotai gives the ratio N signal /Ntotai ■ 
Fig. 8 shows the momentum distributions for experimental data and MC predictions. Figs. 
6 and 8 are used for this estimation. 

The function fo{b) in Eq- HH describes the fractions of D mesons in the c(6) decays, and 
the ratio fc{b)/{fc + fb) gives the probability that D candidate is from a primary c(6) quark. 
We regard fc(b) as a function of soft-pion momentum, P. The function is expressed as 
fc(b) = ^c{b) ■ dc{b){P)- Here dc{b) is determined by the shape of MC soft-pion momentum 
distributions in c{b) -^ D and uj^b) is the estimated total fraction of the c(6) -^ D among 
the selected candidates. 

Performing the maximum likelihood fit to the data sample, we measure A^ = 0.654±0.125 
(1993-95) and Ac = 0.673 ± 0.056 (1996-98). As a check, we also measure Ac with a simple 
binned fit as Ac = 0.520 ± 0.164 (1993-95) and Ac = 0.665 ± 0.085 (1996-98), which are 
consistent with the above values. 

The first- and second-order QCD correction and QED correction are applied with the 
same method as in the exclusive D reconstruction analysis. In the QCD correction, the 
correction factor due to the analysis bias is estimated as Cc = 0.40 ± 0.14 for c-quark and 
Cb = 0.19 ± 0.09 for 6-quark. Applying the first- and second-order QCD correction with this 
factors, and QED correction, we obtain Ac = 0.669 ±0.127 (1993-95) and Ac = 0.689 ±0.057 
(1996-98). 



C. Systematic errors 

The estimated uncertainties in this analysis are summarized in Table |I|, where we show 
average systematic errors for the 1993-98 data. In the soft-pion analysis, we use the same 
procedures to estimate the systematic errors as those in the exclusive D^*^ reconstruction 
analysis in many sources. Here we only explain error sources where we take a different 
method. 

The largest uncertainties are due to the imperfect knowledge of the background fraction 
and its shape. The background is determined by fitting to the P^ distribution of the exper- 
imental data, and we try two functions Fi and F2 described above. In order to estimate the 
background fraction uncertainty, we fix the background shape as Fi, and change its height so 
as to cover the possible range of the background fraction. The background shape uncertainty 
is estimated by using the two background shapes, Fi and F2, while keeping the integrated 
number of the background events in the signal region ( P^ < 0.1 (GeV/c)^) constant. 

The shape of the soft-pion momentum distributions in b ^ D* or c -^ D* is determined 
by fitting to the MC distributions. The uncertainty concerning this distribution is estimated 
by performing the analysis using a binned momentum distribution instead of fitting. 

The total systematic errors are obtained to be ±0.067 and ±0.053 for 1993-95 and 1996- 
98, respectively. 

D. Results 

The Ac values obtained in the inclusive soft-pion analysis are A^. = 0.669 ±0.127(stat.) ± 
0.067{sys.) (1993-95) and A^ = 0.689 ± 0.057(stat.) ± 0.053{sys.) (1996-98). The combined 
result is 

Ac = 0.685 ± 0.052(stat.) ± 0.038{sys.) 



V. CONCLUSION 

Using the 1993-98 experimental data collected by the SLD experiment, we measure the 
parity-violation parameter Ac using two different c-quark tagging methods: 

Ac = 0.690 ± 0.042(stat.) ± 0M9{sys.) and 

Ac = 0.685 ± 0.052(stat.) ± 0.036{sys.), 

from exclusive charmed-meson reconstruction and inclusive soft-pion analysis, respectively. 
To combine them, we must avoid double counting signal events from both samples. We 
find that 1182 events are common to the two analyses. The statistical error for the soft-pion 
analysis without the overlapping events is ±0.061. The combined result is 

A = 0.688 ±0.041, 

where we have also treated the common systematic errors as fully correlated. 

The result is consistent with the standard model prediction of 0.667, obtained by using 
ZFITTER(6.23) with a top-quark mass of 175 GeV/c^, and a Higgs mass of 150 GeV/c^. 
This result represents the currently most precise measurement of Ac- 
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FIG. 1. The distribution of the 2-D impact parameter of the D^ momentum vector to the IP 
for the decay of D*^ -^ D^irf , D^ -^ K^tt^. The sohd circles indicate the experimental data, and 
histograms are MC of D*^ from c-quark (open), from 6-quark (single hatched), and RCBG (double 
hatched). 



SLD 




^cj 140 


— 


^ 


- 


^^ 120 


— 


O 


- 


<N 


_ 


O 100 


— 


o 


_ 


^ 80 


- 


\ 


- 


g 60 


- 






f^ 




e 40 


- 


hq 


- 


20 


- 




Am (GeV/c) 



0.12 0.14 0.16 0.18 0.2 

Am (GeV/c^) 




Am (GeV/c) 



0.14 0.16 0.18 0.2 

Am (GeV/c^) 



FIG. 2. The mass-difference distributions for the decay of {a)D*^ 



D'^TT+, D° -^ K-Tr+, 



(h)D^ -^ K~TT+TT°, (c)Z)0 -^ K-TT+TT+TT-, and (d)Z)0 -^ K'l+i^i (/ =e or fi). The sohd cir- 
cles indicate the experimental data, and histograms are MC of signal (open) and RCBG (double 
hatched). 
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FIG. 3. The mass distributions for {a)D^ and {h)D^ mesons. The sohd circles indicate the 
experimental data, and histograms are the MC of signal (open) and RCBG (double hatched). The 
peaks around m^Kn) '^ 1.6 GeV/c^ in figure (b) comes from the decay D^ -^ Ktttt^. 
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FIG. 4. The distributions of g-cos Oo for the selected D meson sample for (a) left- and (b) right- 
handed electron beams. The solid circles are experimental data, and double hatched histograms 
are RCBG estimated from side-band regions. 
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FIG. 5. The xd distributions for {a)D*^ 
(c)L>o -^ K-TT+TT+TT-, (d)D° ^ R-l+i^i {I =6 or n) {e)D+ -^ K-'k+tt+, and (f)D° -^ K--K+. 
The soUd circles are experimental data and hatched histograms are background estimated from 
side-band events. MC predictions for D mesons from c-decay (open histograms) and 6-decay (sin- 
gle hatched histogram) are also shown. 
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FIG. 6. The P^ distributions for soft-pion candidate tracks, (a) The sohd circles indicate the 
experimental data. The curves are the result of the a fit S{P^) + Fi(Pj,) performed for Pj, < 0.1 
GeV/c (solid line), and the extrapolations of Fi{P^) (dashed line) and F2{P^) (dotted line). The 
definition of the functions are described in the text, (b) The solid circles are the experimental data, 
and histograms are MC predictions for D mesons from c-decay (open), D mesons from 6-decay 
(single hatched), and background (double hatched). The extrapolation of Fi{P^) is also shown as 
a dashed line. 
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FIG. 7. The distributions of q ■ cos^d for the selected D*+ meson sample for (a) left- and (b) 
right- handed electron beams. The solid circles are experimental data, and hatched histograms are 
RCBG estimated from side-band regions. 
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FIG. 8. The momentum distribution for soft-pion candidate tracks. The points are experi- 
mental data. The histograms are MC predictions of D's from c decays (open), D's from b decays 
(single hatched), and background (double-hatched). 



TABLES 
TABLE I. The number of selected candidates from 1993-98 SLD experimental data, and con- 
tributions from c ^ D, b ^ D, and RCBG estimated by MC. 



Channel 


Candidates 


c - 


^ D 


b^ D 


RCBG 


D*+ ^ D^irt, 












DO -^ K-TT+ 


561 


413 


(74%) 


59 (10%) 


89 (16%) 


DO -^ K-tt+ttO 


896 


601 


(67%) 


83 ( 9%) 


212 (24%) 


D^ -^ K-TT+TT+TT- 


537 


418 


(78%) 


36 ( 7%) 


83 (15%) 


D^ ^K-l+9 


433 


296 


(68%) 


31 ( 7%) 


106 (24%) 


D+ -^ K-TT+TT+ 


957 


698 


(73%) 


45 ( 5%) 


214 (22%) 


DO -^ K-TT+ 


583 


403 


(69%) 


27 ( 5%) 


153 (26%) 


Total 


3967 


2829 


(71%) 


281 ( 7%) 


857 (22%) 



TABLE II. 1993-98 average contributions to the estimated systematic error for exclusive D 
meson reconstruction analysis (left column) and inclusive soft-pion analysis (right column). 







SAc 


Source 


Exclusive D^*^ 


Inclusive soft pion 


Background fraction 


0.0111 


0.0324 


Background acceptance 


0.0087 


0.0122 


Background xd / P^ distribution 


0.0112 


0.0018 


Background asymmetry 


0.0028 


0.0093 


fb-*D/ifb~*D + Ic^d) 


0.0011 


0.0018 


Ab-,D (Ah) 


0.0017 


0.0021 


Ab^D (Mixing) 


0.0092 


0.0120 


c fragmentation 


0.0003 


0.0010 


b fragmentation 


0.0003 


0.0005 


D meson x^ shape / Soft-pion momentum shape 


0.0040 


0.0003 


Polarization 


0.0035 


0.0033 


Ae 


0.0002 


0.0005 


as 


0.0004 


0.0005 


Correction factor for first order QCD correction 


0.0024 


0.0033 


Second order QCD correction 


0.0006 


0.0008 


Gluon splitting 


0.0002 


0.0005 


Total 


0.0213 


0.0383 



